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Preface 

) 

This  thesis  grew  out  of  a study  done  by  AFIT  faculty 
members  Stan  Robinson  ajid  Jurgen  Gobien  which  proposed 
statistical  model  for  the  infrared  background  field  and  -^he 
resulting  detector  output.  The  receiver  structures  developed 
in  this  thesis  are  based  upon  this  model.  There  are  many 
complex  problems  associated  with  the  performance  problem 
addressed  by  this  study  and  time  constaints  prevented  me 
from  examining  them  all.  I hope  this  study  will  at  least 
provide  some  insight  and  possibly  serve  as  a base  from 
which  further  studies  can  develop. 

I would  like  to  express  my  appreciation  to  my  thesis 
advisor  and  instructor,  Stan  Robinson,  who  taught  me  more 
I about  random  processes  and  statistics  than  I really  wanted 

to  know.  Although  I've  learned  alot  in  my  studies  here  at 
AFIT,  at  this  point  I am  glad  that  it  is  now  in  the  past. 

Stephen  J.  Dunning 
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Abstract 

The  use  of  passive  electro-optical  sensors  in  military- 
systems  has  received  increased  emphasis  in  recent  years.  In 
particular,  passive  infrared  sensors  could  potentially  be 
used  in  an  infrared  airborne  threat  warning  receiver  which 
would  detect  the  infrared  signature  of  a hostile  threat  and 
provide  a warning  indication  to  the  aircraft  crew. 

The  performance  of  an  optical  receiver  designed  to  de- 
tect a target  by  taking  advantage  of  the  target's  spectral 
signature  is  presented.  The  receiver  processes  the  signal 
in  several  narrow  frequency  bands  and  is  based  upon  a sta- 
tistical model  which  represents  the  field  in  each  band  as  a 
Gaussian  random  process  whose  moments  depend  upon  the  target 
and  background  characteristics.  The  signal  is  detected  by 
am  array  of  power  detectors  whose  outputs  are  modeled  as 
random  variables  characterized  by  non-central  chi  square 
probability  density  fimctions. 

The  optimal  Bayes/Neyman-Pearson  receiver  structure  for 
an  M spectral  channel,  N sequential  look  scanning  receiver  is 
presented  and  two  practical  suboptimal  receiver  structures 
are  developed.  An  attempt  is  made  to  obtain  closed  form  an- 
alytical performance  expressions,  but  it  quickly  becomes 
evident  that  closed  form  expressions  are  obtainable  only  for 
a few  special  cases.  Thus  numerical  methods  using  a digital 
computer  are  used  to  calculate  the  probability  of  false  alarm 
and  the  probability  of  detection  for  each  processor  using 
identical  parameters.  These  results  are  used  to  plot  a re- 


ceiver  operating  curve  (ROC)  for  each  processor  and  are  com- 
pared to  the  performance  of  two  ad-hoc  linear  processors# 
ROC'S  are  obtained  for  variations  in  parameters  (mean  and 
variance)  and  for  one,  two,  and  three  spectral  channel 
receivers . 

The  results  indicate  that  the  approximate  receiver 
structures  and  the  ad-hoc  receiver  structures  all  have  the 
same  performance.  This  is  demonstrated  analytically  for  the 
case  where  only  the  mean  is  varied.  The  results  also  show 
that  performance  depends  only  upon  the  difference  in  the 
square  roots  of  the  mean  to  variance  ratios  under  each 
hypothesis  and  the  ratio  of  the  variances.  It  is  concluded 
by  this  study  that  performance  equal  to  that  of  either 
receiver  structure  derived  from  the  general  optimal  processor 
can  be  obtained  by  using  a relatively  simple  ad-hoc  linear 
receiver. 
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MULTICHANNEL  INFRARED  RECEIVER 
PERFORMANCE 

I . Introduction 

The  use  of  passive  electro-optical  sensors  has  received 
increased  emphasis  by  the  Air  Force  in  recent  years.  Pas- 
sive optical  and  infrared  sensors  have  been  developed  for  a 
nvuaber  of  military  applications.  In  particular,  passive 
infrared  sensors  could  potentially  be  used  in  an  infrared 
airborne  threat  warning  receiver  which  would  detect  the 
infrared  signature  of  a hostile  threat  and  provide  a warning 
indication  to  the  aircraft  crew. 

Background 

The  use  of  airborne  threat  warning  receivers  employing 
scanned  infrared  detectors  has  often  been  proposed  as  a 
method  of  detecting  airborne  threats  such  as  anti-aircraft 
missiles.  The  practicality  of  such  scanning  systems  is 
supported  by  the  well  developed  sensor  technology  aquired 
during  the  Forward  Looking  Infrared  Receiver  (FLIR)  system 
program.  A warning  receiver  on  board  an  aircraft  ideally 
would  detect  a launched  airborne  threat  and  cause  some  form 
of  defensive  countenneasure  to  occur. 

As  avionics  systems  become  more  automatic  in  their 
response  to  the  aircraft  environment,  there  is  a obvious 
need  to  reduce  the  false  alarm  rate  of  any  eventual  airborne 
threat  warning  receiver  while  maintaining  acceptable  threat 
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detection  performance.  A large  number  of  false  alarms  would 
generate  di strust  in  the  human  operator  and  would  also  un- 
necessarily expend  countermeasure  resources.  A well  designed 
threat  warning  receiver  would  have  to  use  signal  processing 
techniques  which  would  enable  the  receiver  to  distinguish 
between  valid  targets  and  phenomenon  such  as  sun  glint,  open 
fires,  smoke  stacks,  and  other  non-hostile  thermal  sources 
which  could  contribute  to  false  alarms. 

Problem  Statement 

The  purpose  of  this  thesis  is  to  apply  statistical 
signal  detection  techniques  to  a statistical  infrared  back- 
groimd  model  proposed  by  Gobien  and  Robinson  (Ref  11)  to 
determine  the  structure  and  performance  of  an  optimal  multi- 
channel infrared  receiver. 

Approach  and  Scope 

The  receiver  structure  and  performance  problem  addressed 
by  this  thesis  was  divided  into  three  distinct  areas;  the 
statistical  model,  the  receiver  structure,  and  receiver 
(signal  processor)  performance. 

Statistical  Model.  Any  receiver  structure  obtained 
through  the  use  of  statistical  signal  detection  techniques 
requires  a statistical  model  which  describes  the  signals  the 
receiver  is  to  detect.  For  the  purposes  of  this  study,  the 
signals  to  be  processed  by  the  receiver  were  the  presence 
of  a valid  target  (the  target  or  one  hypothesis)  and  the 
absence  of  a valid  target  (the  null  or  zero  hypothesis). 
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The  statistical  model  proposed  by  Gobien  and  Robinson  and 
used  throughout  this  thesis  states  that  the  infrared  back- 
groxmd  field  may  be  represented  by  a Gaussian  random  process 
whose  moments  depend  upon  the  target  and  background  charac- 
teristics. The  model  further  states  that  the  output  of  a 
power  detector  that  intercepts  this  field  is  characterized 
by  a non-central  chi-square  probability  density  function. 

Receiver  Structure.  Once  the  statistical  signal  model 
was  assumed,  statistical  signal  detection  techniques  were 
applied  to  the  model  to  arrive  at  a binary  receiver  structure 
capable  of  distinguishing  between  the  two  previously  defined 
hypotheses.  The  likelihood  ratio  was  formed  and  an  optimal 
Bayes/Neyman-Pearson  receiver  processor  structure  was 
obtained.  Because  of  the  complex  structure  of  the  optimal 
processor,  large  argument  and  small  argument  approximations 
were  used  to  simplify  the  signal  processor  structure.  These 
approximations  resulted  in  two  practical  processor  structures, 
a non-linear  receiver  processor  and  a small  signal  linear 
receiver  processor. 

Receiver  Performance.  To  determine  receiver  performance, 
an  attempt  was  first  made  to  determine  analytic  performance 
expressions  such  as  the  Chemoff  bound  for  each  processor 
structure.  It  quickly  became  evident  that  the  complexity  of 
the  mathematics  involved  would  prevent  the  derivation  of  any 
useful  analytic  results  and  necessitated  the  use  of  numerical 
methods  to  determine  performance.  Receiver  performance  was 
calculated  numerically  using  the  Control  Data  CYBER-7^ 
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computer  system  available  to  AFIT.  FORTRAN  programs  were 
written  to  compute  the  processor  output  probability  density 
fimctions  for  each  receiver  structure  which  were  then 
numerically  integrated  to  yield  receiver  performance  in 
terms  of  a receiver  operating  curve  (ROC). 

Numerical  analysis  of  receiver  performance  was  limited 
by  time  constaints  to  single  channel  receiver  performance 
for  five  sets  of  parameters,  and  to  two  and  three  chaumel 
receiver  performance  for  a single  set  of  parameters. 
Identical  parameters  were  used  in  the  evaluation  of  each 
receiver  structure,  and  while  the  chosen  parameters  were  not 
necessarily  representative  of  those  one  might  obtain  from  a 
real  system,  they  did  provide  a common  basis  for  comparison 
of  receiver  processor  performance. 
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II.  Statistical  Field  and  Detector  Model 


The  binary  signal  detection  receivers  presented  in  this 
thesis  were  designed  to  detect  the  presence  of  a v^ilid  target 
against  infrared  background  clutter  by  processing  the  output 
signal  from  an  optical  detector.  Before  the  signal  processor 
structure  could  be  determined,  the  statistics  of  the  signals 
to  be  detected  had  to  be  ascertained.  This  chapter  presents 
the  statistical  model  for  the  detector  output  signal  used 
throughout  this  study.  The  statistics  of  the  observed 
optical  field  were  first  examined,  and  these  statistics  were 
then  used  to  derive  the  statistics  for  the  output  signal  of 
an  ideal  power  detector. 

Field  Representation 

The  incident  field  to  be  detected  at  a position  r can 
be  described  by  the  scalar  field  quantity  u(r,t),  which 
may  represent  either  the  electric  or  magnetic  field.  For 
convenience,  this  field  is  represented  by  its  complex 
envelope 

U(r,t)  = UR(r,t)  + jUi(r,t)  (1) 

where  U{^(r,t)  is  the  real  part  of  U(r,t)  and  Uj(r,t)  is  the 
imaginary  part  of  U(r,t).  This  representation  is  the  same 
as  the  quadrature  field  model  used  in  radar  and  communication 
systems.  The  complex  envelope  is  implicitly  defined  by  the 
relationship 

u(r,t)  = Re[u(r,t)  exp(- j2nfot)]  (2) 
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where  Tq  is  the  optical  center  frequency  and  Re(*)  denotes 
the  real  part  of  the  enclosed  quantity  (Ref  35I8IO).  Because 
the  complex  envelope  U(r,t)  is  a time  varying  quantity,  it 
is  a valid  representation  for  extremes  ranging  from  the 
incoherent  light  due  to  naturally  occurring  illumination,  to 
the  coherent  output  of  a laser.  For  the  purpose  of  this 
thesis,  the  complex  envelope  was  used  to  represent  the 
received  fields  in  the  intermediate  infrared  region. 

It  is  unreasonable  to  assume  a priori  the  field  that  the 
detector  would  intercept.  Therefore,  it  is  appropriate  to 
think  of  the  complex  envelope  of  the  field  as  a random  pro- 
cess in  both  time  and  space.  It  was  further  assumed  that 
the  complex  envelope  was  a sample  fimction  of  a complex 
Gaussian  raindom  process.  This  assumption  can  be  justified 
by  the  application  of  the  central  limit  theorem  (Ref  10:266) 
to  the  received  field,  where  the  received  field  is  due  to 
the  sum  of  a large  number  of  individual  fields,  each  of 
which  is  due  to  the  scattering  of  natural  light  by  an 
independent  particle . 

For  simplicity,  the  real  and  imaginary  parts  of  U(r,t) 
are  assumed  to  be  independent  and  identically  distributed 
(Ref  3:1849).  The  complete  specification  is  then  given  by 

E[u(r,t)]  = m(r,t)  (3) 

E[U(r,t)  U«(r',t')]=  R(r,r',t,f)  (4) 

where  the  above  conditions  imply  that  R(r,r',t,t')  is  a real 
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ftinction  and  is  twice  the  correlation  function  of  either  the 
real  or  imaginary  part  of  the  field. 

The  optical  detection  problem  can  now  be  stated.  When 
a target  is  present,  the  received  field  becomes 

U(r,t)  = Us(r,t)  + Ub(r,t)  (5) 

where  Ug(r,t)  is  the  signal  or  target  field  and  U^(r,t)  is 
the  background  field.  This  hypothesis  is  denoted  by  . 

If  no  target  is  present,  denoted  by  Hg  i the  received  field 
is  only  U^(r,t).  Under  the  assumption  that  Ug(r,t)  and 
Ui3(r,t)  are  independent,  the  statistics  under  each  hypothesis 
are  given  by 

m(r,t)  = mt)(r,t) 

Ho* 

R(r,r',t,t')  = Rij(r,r' ,t,t' ) (6) 


m(r,t)  = ms(r,t)  + mt,(r,t) 

R(r,r',t,t')  = Rs(r,r’ ,t,t' ) + Rt)(r,r' ,t,  t' ) (7) 

With  the  field  model  described  above,  it  is  possible  to 
completely  describe  the  space-time  processing  that  should  be 
accomplished  by  a receiver  and  the  resulting  performance 
(Ref  13).  There  are  three  major  disadvantages  to  the  space- 
time  processing  approach  (Ref  11).  The  processing  and  per- 
formance depend  upon  m(r,t)  and  R(r,r',t,t')  explicitly,  and 
it  is  unreasonable  to  expect  that  these  quantities  are  known. 
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Secondly,  the  space- time  processing  required  would  typically 
be  much  too  complex  to  implement  even  if  m(r,t)  and  R(r,r ' , 
t,t')  were  known.  Finally,  the  processing  would  require  the 
measurement  of  the  complex  envelope  U(r,t)  and  this  is  not 
possible  at  most  wavelengths  using  currently  available 
devices.  Although  equations  (6)  and  (?)  completely  specify 
the  field  statistics,  the  determination  of  m(r,t)  and 
R(r,r',t,t')  explicitly  is  generally  quite  difficult.  Thus, 
reasonable  approximations  were  sought  which  would  permit 
more  practical  and  realizable  signal  processor  configurations. 

Two  simple  but  crude  parameters  which  were  used  to 
simplify  the  field  statistics  were  coherence  distance  and 
coherence  time,  denoted  D^,  and  T^  , respectively.  The 
coherence  distance  is  described  by 

R(r,r',t,t  ) = m(r,t)  m(r',t)  ; |r  - r'|>Dc  (8) 

for  all  t.  It  is  the  distance  beyond  which  samples  of  the 
field  are  considered  to  be  statistically  independent.  The 
resulting  coherence  cell  model  for  the  field  is  a simple 
approximation  in  which  it  is  assumed  that  the  incident  field 
is  spatially  constant  over  an  area  (coherence  cell)  and  is 
statistically  independent  from  the  field  in  other  cells. 

Coherence  time  is  defined  as  the  length  of  time  over 
which  the  field  can  be  broken  into  piecewise  constant  inde- 
pendent samples  in  time  and  is  described  by 
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(9) 


R(r,r, t, t+T(,)  = m(r,t)  m(r,t+T(j) 

for  all  r.  The  coherence  time/cell  model  allows  the  field 
to  he  decimated  in  space  and  time  so  that  the  field  cam  now 
he  considered  as  a collection  of  random  variables  rather 
tham  a ramdom  process  in  space  and  time. 

The  field  can  also  he  separated  into  a number  of  disjoint 
frequency  windows.  Since  the  fields  are  considered  to  he 
uncorrelated,  the  output  field  in  any  spectral  window  is 
independent  from  the  field  in  any  other  frequency  window. 

This  independence  between  frequency  windows  allows  the 
coherence  time/cell  field  model  developed  above  to  he  applied 
to  the  statistical  description  of  the  field  in  each  spectral 
window,  differing  only  in  the  moments  required  to  describe 
each  output  field. 

Detector  Model 

Utilizing  the  piecewise  constamt  coherence  time/cell 
model  for  each  spectral  window,  an  array  of  detectors  placed 
in  the  measurement  plane  of  the  receiver  can  he  considered. 
Each  detector  is  an  ideal  power  detector  whose  active  area 
is  matched  to  the  smallest  coherence  cell  expected,  where 
the  size  of  the  coherence  cell  depends  upon  the  type  of 
target  to  be  detected.  Coherence  cell  size  in  excess  of  that 
required  would  lead  to  performance  degradation  due  to  the 
increase  in  background  noise. 

The  intensity  or  rate  function,  A.(t),  of  the  j"^^  detector 
is  given  by 
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A.4(t)  = ^ / |u(r,t)|  ^ dr 


(10) 


where  n is  the  quantum  efficiency  of  the  detector,  hfo  is  the 
energy  of  a photon,  is  the  active  area  of  the  detector, 

and  ■ t is  assumed  that  the  intensity  of  the  scalar  field  in 

units  of  power  per  unit  area  is  given  by  |u(r,t)|  The  rate 
f\mction  may  be  defined  as  the  average  number  (ensemble)  of 
photo-electrons  observed  at  the  output  of  the  detector  as  a 
function  of  time. 

Within  the  constraints  of  the  coherence  time/cell  model 
and  assuming  that  the  detector  area  Aj  is  on  the  order  of  the 
coherence  cell  and  that  the  time  interval  in  which  the  obser- 
vation is  made  is  less  than  the  coherence  time,  the  current 

output  of  the  detector,  centered  at  rj  , is  given  by 

y.  = qX.^(t)  = ^ J'|U(r,t)|  ^ dr 
J J hlo 

where  q is  the  charge  of  an  electron. 

The  statistics  of  the  detector  output  can  now  be  deter- 
mined. The  randomness  of  the  detector  output  y-  is  due  only 

J 

to  the  stochastic  nature  of  the  received  fields,  since  the 
ideal  nature  of  the  detector  eliminates  any  device  noise 
that  would  be  inherent  in  a real  photodetector.  The  ideal 
nature  of  the  detector  serves  to  simplify  the  analysis  and  is 
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based  upon  the  contention  that  detector  noise  need  not  be 
considered  if  the  noiseless  performance  of  the  receiver 
proved  to  be  unacceptable.  The  probability  density  function 
(pdf)  of  the  output  current  of  the  detector  is  given  by 


(12) 


where  IqC*)  is  the  modified  Bessel  function  of  the  first  kind 
of  order  zero.  This  pdf  is  know  as  the  non-central  chi- 
square  density  function  with  two  degrees  of  freedom.  The 
mean  and  variance  of  this  pdf  are  related  to  the  field  quan- 
tities by 

(13) 

(14) 


'j  = 2 Aj  m(rj,t) 


2 _ / qn 


hfo 


2 

- (m(rj,t)) 


It  is  straightforward  to  extend  Eq  (12)  to  the  joint 
probability  density  function  for  the  detector  array.  Since 
the  detectors  are  disjoint  in  frequency,  and  therefore 
independent,  the  joint  pdf  is  given  by 

M 

f(l)  = n fj-(yj-)  (15) 
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Eqs  (12)  and  (15)  also  apply  to  the  detector  output  of  a 
single  detector  which  is  scanned  across  the  image  plane  a 
coherence  cell  at  a time.  If  the  frame  time  for  each  coher- 
ence cell  is  on  the  order  of  the  coherence  time,  the  detector 
output  for  each  frame  time  is  independent  and  the  joint  pdf 
is  given  by  Eq  (15) i where  the  index  j refers  to  the  time 
frame  (or  coherence  cell)  in  which  the  measurement  is  made. 

Infrared  Background  Data 

To  investigate  the  validity  of  the  detector  output  model 
developed  above,  a literature  search  was  conducted  to  dis- 
cover any  sources  of  experimental  data  pertaining  to  the 
statistical  properties  of  infrared  backgrounds.  While  many 
sources  relating  the  spectral  characteristics  of  infrared 
backgrounds  were  found,  only  two  sources  concerning  the 
statistical  properties  of  infrared  backgrounds  were  discovered. 

The  first  source  of  statistical  information  was  found 
in  a 197^  paper  titled  "Statistical  Properties  of  the  Back- 
groxind  Noise  for  the  Atmospheric  Windows  in  the  Intermediate 
Infrared  Region"  (Ref  4).  Measurements  were  made  of  various 
background  types  using  a scanning  radiometer.  The  authors 
analyzed  the  data  and  developed  a statistical  model  for  the 
infrared  background.  While  the  model  proposed  in  the  paper 
is  based  upon  a combination  of  Gaussian  and  Poisson  statistics 
and  is  different  from  the  model  used  in  this  thesis,  the 
statistical  data  presented  does  support  the  non-central  chi- 
square  model.  The  distinctive  shapes  of  the  probability 
density  functions  illustrated  in  the  report  (Ref  4j28)  are 


12 


all  possible  forms  of  the  non-central  chi  square  probability 
density  function,  dependent  upon  the  specific  parameters  of 
the  density. 

Another  source  of  infrared  background  statistical  data 
was  found  in  a series  of  reports  by  the  Lockheed  Missile  and 
Space  Company  (LMSC)  (Refs  5;  6;  7;  and  8).  Under  the  spon- 
sorship of  the  Advanced  Research  Projects  Agency  and  the  U.S. 

Army  Missile  Command,  LMSC  began  the  Backgroxmd  Measurements 
Program  in  which  natural  infrared  backgrounds  were  measured 
from  the  air  using  an  infrared  radiometer  mounted  in  a U-2 
research  aircraft. 

In  one  report,  data  for  one  of  the  measurement  flights 
was  analyzed  by  LMSC  and  the  data  presented  in  the  form  of 
histograms  for  each  of  six  spectral  filters  (Ref  5)»  Tbe 
characteristics  of  the  six  filters  are  listed  in  Table  I. 

While  the  histograms  contain  data  for  a combination  of 
different  background  types  that  were  overflown  by  the  air- 
craft, comparison  of  the  data  with  the  flight  track  by  LMSC 
indicated  that  it  was  possible  to  separate  the  histograms  for 
each  of  four  background  types  from  the  combined  histogram. 

The  four  background  types  were  high  cloud,  low  cloud,  water, 
and  terrain.  The  experimental  histogram  for  filter  1 is 
shown  in  Figure  1 . 

If  the  number  of  sample  points  is  large,  the  histogram 
or  sample  relative  frequency  plot  for  each  background  type 
will  converge  to  the  corresponding  pdf  for  that  background. 

The  weighted  sum  of  these  pdf's  will  yield  a cumulative  rela- 
tive frequency  plot  or  histogram.  To  test  the  fit  of  the 

I 
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Table  I 

Filter  Characteristics 


Filter  No. 

Center  Wavelength  (pm) 

Bandwidth  (pm) 

1 

4.50 

0.202 

6 

4.48 

0.157 

5 

4.4^ 

0.124 

4 

4.44 

0.097 

3 

4.42 

0.074 

2 

4.4l 

0-055 

(From  Ref 

model  with  the  experimental  data,  non-central  chi-square 
probability  density  functions  were  generated  for  each  back- 
ground type  using  parameters  estimated  from  the  histogram  in 
the  LMSC  report.  The  pdf's  generated  for  filter  1 for  each 
background  type  are  illustrated  in  Figures  2 - 5»  These 
pdf's  were  then  linearly  combined  according  to  the  following 
equation  which  was  derived  through  trial  and  error: 

f(x)  = 0.08  f}{Q(x)  + 0.14  fj^(x) 

+ 0.26  f^(x)  + 0.52  f^(x)  (16) 

where  f(x)  is  the  combined  relative  frequency  function,  fjjQ(x) 
is  the  high  cloud  background  pdf,  f]jQ(x)  is  the  low  cloud 
background  pdf,  fyj(x)  is  the  water  background  pdf,  and  fqi(x) 
is  the  terrain  background  pdf.  The  numerical  constants  in 
Eq  (16)  are  the  estimated  fraction  of  total  samples  con- 
tributed to  the  combined  relative  frequency  function  by  each 
type  of  background.  The  combined  relative  frequency  plot 
obtained  for  filter  1 is  illustrated  in  Figure  6.  This  sample 
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RADIANCE  IN  MICR0WATTS/CCM2-SR-MICR0ND 


Fig.  1.  Experimental  Histogram,  Filter  1 (From  Ref  5:4-2) 


4 


relative  frequency  plot  closely  approximates  the  experimental 
histogram  obtained  by  LMSC . 

This  procedure  was  repeated  for  filters  4,  5»  £Ln<i  6. 
Density  parameters  were  estimated  for  each  background  type 
for  each  filter  and  the  pdf's  which  were  obtained  were  com- 
bined according  to  Eq  (l6)  to  obtain  a cumulative  relative 
frequency  plot  for  each  filter.  These  plots  are  illustrated 
in  Appendix  A and  also  closely  approximate  the  histograms 
obtained  experimentally  by  LMSC  for  the  corresponding  filter. 
The  data  for  filters  2 and  3 was  so  uniform  that  reliable 
estimates  of  pdf  parameters  could  not  be  made. 

While  the  experimental  data  reviewed  here  does  not  con- 
clusively prove  the  validity  of  the  non-central  chi-square 
model,  it  does  indicate  that  the  model  is  consistent  with 
the  experimental  infrared  background  statistical  data 
currently  available. 
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NON-CENTRRL  CHI  SQ  DENSITY 


Fig.  2.  Estimated  pdf,  High  Cloud  Background 
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NON-CENTRflL  CHI  SQ  DENSITY 


Fig.  3*  Estimated  pdf,  Low  Cloud  Background 


18 


NON-CENTRflL  CHI  SO  DENSITY 


Fig.  4.  Estimated  pdf,  Water  Background 
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NON-CENTRflL  CHI  SQ  DENSITY 


Fig.  5»  Estimated  pdf,  Terrain  Background 
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III.  Signal  Processor  Structures 

A signal  processor  structure  developed  through  the  appli- 
cation of  statistical  signal  detection  techniques  depends 
explicitly  upon  the  statistical  characteristics  of  the  sig- 
nals being  detected.  The  signals  to  be  processed  by  the  pro- 
cessors developed  here  are  the  detector  outputs  under  the 
null  and  target  hypotheses.  By  correctly  processing  these 
signals,  it  is  hoped  that  the  receiver  will  distinguish 
between  the  presence  or  absence  of  a valid  target  with  a high 
degree  of  accuracy.  This  chapter  presents  the  receiver  pro- 
cessor structures  developed  by  using  the  non-central  chi- 
square  detector  output  model  to  characterize  the  statistics 
of  the  null  and  target  hypothesis  signals. 

The  basic  receiver  considered  here  is  a scanning  receiver 
which  , by  means  of  narrowband  filters  and  parallel  ideal 
power  detectors,  observes  M frequency  disjoint  channels. 

The  reasoning  behind  this  structure  lies  in  the  fact  that 
most  targets  to  be  detected  will  have  characteristic  spectral 
signatures.  By  properly  choosing  the  spectral  channels,  the 
receiver  can  discriminate  against  unwanted  sources  whose 
spectral  characteristics  differ  from  those  of  the  desired 
target. 

The  receiver  processor  development  and  following  analy- 
sis pertains  to  a receiver  which  makes  N successive  obser- 
vations of  M detector  outputs  for  an  arbitrary  coherence  cell. 
The  channel  outputs  during  the  i'^h  observation  of  the  cell 
would  constitute  a vector  = Yu  »yi2' • • • *yiM  where  the 
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elements  are  the  chainnel  output  currents.  Each  successive 
observation  is  assumed  statistically  independent  and  iden- 
tically distributed.  The  first  assumption  is  based  upon  the 
coherence  time/cell  model,  while  the  second  assumption  serves 
to  simplify  the  processor  structure  aind  later  analysis  by 
excluding  temporal  processing  of  successive  observations. 

If  the  target  has  known  temporal  characteristics,  temporal 
processing  would  be  advantageous,  but  at  the  cost  of  increased 
receiver  complexity. 

Optimal  Processor  Structure 

The  totality  of  all  observations  made  by  the  receiver  is 
denoted  by  y and  has  the  following  joint  probability  density 
functions  under  the  null  and  target  hypotheses: 


^ N M . 

n n — 7-  exp 
k=l  i=l  2CTio 


Ho‘ 


X I, 


72“ 

C^iO 


ki  ™i0 


^ki 


4-0 


0 ; elsewhere 
N M 1 

n n — 2~ 

k=l  i=l 


(17) 


^ (Yki+niil) 
2^11 


fi(y)  =< 


X I, 


— 2“  V^ki  ”“il' 
Lc^il 


! yki-0 


0 j elsewhere 


(18) 
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These  pdf's  are  obtained  from  Eq.s  (12)  and  (15)'  The  receiver 
must  choose  between  the  hypotheses  by  processing  the  received 
signal  described  by  the  above  probability  density  functions. 

An  optimal  decision  is  made  by  comparing  the  likelihood  ratio, 
A(y)  = (y)/f‘o(y) » with  a threshold  and  declaring  a target 

if  the  threshold  is  exceeded  (Ref  12:19-^6).  An  equivalent 
expression  is  obtained  by  computing  the  logarithm  of  the  like- 
lihood ratio.  This  results  in  the  statistic  z(y),  which  is 
defined  by 

N M 

z(y)  = In  A(y)  =22  Yki^  lo^^iV^) 

k=l  i=l ^ 

^ - In  Io(®iVyki'^J  (19) 

where 

2 2 

3.^  = 2 ln(crj^Q/CTj^l ) + i[(niio/'^iO^  " ^’^il/'^il^] 

bi  = i(l/afo  - Vctii) 

Ci  = 

®i  = V^SIo/'^iO 

If  the  threshold  is  v,  then  the  optimal  signal  processor 
algorithm  is  given  by 

N M 1^1 

2 2 [Di  In  - In  I„(eii,'yJJ)]  $ w (20) 

K-1  1-1  Dq 
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M 

where  w = ln(v)  - NE 

i=l  ^ 

and  and  Dq  represent  the  decisions  "a  valid  target  is 
present"  and  "no  target  is  present,"  respectively.-  One 
cheinnel  of  the  required  processor  structure  is  illustrated 
in  Figure  ?•  The  complex  processing  required  makes  this 
structure  impractical  for  a real  time  system  and  also  makes 
performance  analysis  extreniely  difficult. 


Fig.  7*  Optimal  Signal  Processor  Structure 


25 


Non-Linear  Approximate  Processor 

The  optimal  processor  structure  may  be  simplified  by 
substituting  the  large  argximent  approximation  for  the  modi- 
fied Bessel  fimction  into  the  optimal  signal  processor 
algorithm  given  in  Eq  (20).  This  approximation  is  given  by 

Iq(x)  « e^/^f2nx  ; x>l  (21) 

Substituting  this  approximation  into  the  processor  algorithm 
and  reducing  the  expression,  the  new  signal  processor  algo- 
rithm becomes 


N 

2 

k=l 


M 

.2  yki+(Ci-  ei)Vj^ 

1=1  Dq 


M 

where  w = ln(v)  - N 2 i ln(ei/ci ) 

i=l  ^ 


(22) 


The  processing  required  by  this  structure  is  less  complex 
than  that  required  by  the  optimal  processor  and  is  much  more 
practical  to  implement.  One  channel  of  the  non-linear  pro- 
cessor structure  is  illustrated  in  Figure  8. 

Linear  Approximate  Processor 

By  approximating  the  modified  Bessel  function  with  its 
small  argument  equivalent,  the  signal  processor  structure 
can  be  further  simplified.  The  small  argument  approximation 
for  the  modified  Bessel  function  is  given  by 
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Fig.  8.  Non-Linear  Approximate  Processor  Structure 


Iq(x)  » xv^  ! X<1  (23) 

Substituting  this  expression  into  the  optimal  processor  algo- 
rithm and  reducing  terms,  the  processor  algorithm  becomes 


N M 


2 y.i  I 

=1  1=1  Dr> 


where  w = 


ln(v)  - N 2 

i=l  ^ 


This  linear  signal  processor  structure  is  simpler  than  either 
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the  optimal  processor  or  the  non-linear  approximate  processor 
and  is  one  of  the  most  elementary  signal  processor  structures 
possible.  One  channel  of  the  necessary  processing  is  illus- 
trated in  Figure  9« 

Ad-hoc  Linear  Processors 

An  ad-hoc  signal  processor  is  one  which  is  obtained 
through  intuition  rather  than  analytical  procedures.  Two 
ad-hoc  linear  signal  processors  are  presented  here  and  will 
be  used  for  later  performance  comparison  with  the  non-linear 
approximate  processor  and  the  linear  approximate  processor. 


Fig.  9*  Linear  Approximate  Processor  Structure 
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The  processor  constants  used  here  are  proportional  to  the 
signal  mean  and  inverse  to  signal  fluctuation.  The  algo- 
rithms for  the  two  ad-hoc  processors  are  given  by 
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k=l 
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i=l 


(mil/ail 


mio/a 


^ki 
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"io/^io)  ^ki  ? w 


(25) 


(26) 


where  w is  an  arbitrary  threshold.  The  processor  described 
by  Eq  (25)  will  be  referred  to  as  the  ad-hoc  linear  processor 
#1  and  the  other  processor  will  be  designated  as  the  ad-hoc 
linear  processor  #2.  One  channel  of  the  processor  structure 
for  each  ad-hoc  processor  is  illustrated  in  Figures  10  and  11. 
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IV . Signal  Processor  Performance 


Signal  processor  performance  is  characterized  hy  the 
probability  of  the  receiver  making  an  error.  For  the  Neyman- 
Pearson  processor  structure,  performance  is  completely  spec- 
cified  by  the  quantities  probability  of  detection  (Pp)  and 
probability  of  false  alarm  (PpA^ • The  probability  of 
detection  and  the  probability  of  false  alarm  are  defined  by 


Pe  = P(Di|Hi)  = P[z(y)  > wlHi]  (2?) 

PfA=  P(DilHo)  = P[z(y)  > w|h^  (28) 


These  quantities  are  functions  of  the  threshold  w and  are 
normally  plotted  as  a performance  curve,  Pp  vs  Pfa»  knovm  as 
the  receiver  operating  characteristic  or  receiver  operating 
curve  (ROC).  The  exact  computation  of  the  probability  of 
detection  and  the  probability  of  false  alarm  requires  that 
the  probability  density  function  (or  distibution  function) 
for  the  log  likelihood  ratio,  2(y),  be  knovm  under  both 
hypotheses . 

This  chapter  presents  the  performance  characteristics  of 
the  non-linear  approximate  processor  and  the  linear  approx- 
imate processor  only,  as  the  complexity  of  the  optimal  signal 
processor  makes  the  computation  of  its  output  probability 
density  function  unrealistic.  An  attempt  is  made  to  derive 
analytic  performaince  expressions  by  first  determining  the 
joint  pdf  for  the  output  of  each  signal  processor.  This 
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impossible  except  for  a special  case  of  the  linear  approximate 
processor.  A general  expression  for  the  Chemoff  bound  on 
Pjj  and  Pp^  for  the  linear  approximate  processor  is  derived, 
but  other  analytic  expressions  could  not  be  obtained. 

The  difficulty  in  obtaining  analytic  performance  expres- 
sions for  either  the  non-linear  or  linear  processor  neces- 
sitated the  use  of  niamerical  methods  to  determine  processor 
performance.  Using  a digital  computer,  conditional  processor 
output  pdf's  were  calculated  and  receiver  operating  curves 
for  each  processor  structure  were  generated.  While  numerical 
methods  cannot  yield  absolute  measures  of  performance  unless 
the  exact  "real  world"  density  parameters  are  known,  by 
using  identical  parameters  in  the  performance  calculations 
for  each  processor,  the  relative  performance  of  each  proc- 
essor can  be  determined.  Performance  curves  are  generated 
for  changes  in  mean,  changes  in  variance,  and  for  one,  two, 
and  three  channel  receivers. 

Non-Linear  Approximate  Processor 

Starting  with  the  non-linear  processor  algorithm  given 
by  Eq  (22),  the  processor  output  function,  z(y) , for  a single 
channel  "i"  is  given  by 

z(y)  = bi  y + (ci-ei)Vy’  (29) 

Since  the  detector  output  y is  modeled  by  a non-central  chi- 
square  random  variable,  through  transformation  of  variables 
the  processor  output  probability  density  function  can  be 
derived.  For  any  given  values  of  bf,  cf,  and  ei  there  is  a 
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one-to-one  correspondence  between  y and  z(y),  where  y^O. 
Completing  the  square  and  solving  for  y in  terms  of  z yields 


y(z) 


( V(Ci-ei)  V'bi  + z'  - (Ci-ei)/bi)^ 


(30) 


from  which  the  Jacobian  of  the  transformation  is  found.  The 
Jacobian  is  given  by 


iy(z) 

1 ( ^i”®i ) 

^z 

bi  biV(ci-ei)^  + biz' 

(31) 


where  | • | indicates  the  magnitude  of  the  expression.  The 
unconditional  single  channel  processor  output  probability 
density  function  can  now  be  expressed  by 


f2(z)  = Jfy(  y(z)  ) 
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(32) 


2 

where  mj^  and  ai  are  hypothesis  dependent  channel  parameters 
and  the  range  of  z depends  upon  the  constants  bj^,  Cj^,  and  ej^. 
Since  each  channel  of  a multichannel  receiver  is  disjoint  in 
frequency  and  independent,  the  joint  processor  output  pdf  can 
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be  derived  from  the  product  of  the  characteristic  functions 
for  each  channel. 

The  characteristic  function,  ^(u),  is  a special  case  of 
the  Fourier  transform  given  by  the  relationship 

$(u)  = F(-u/2tt)  (33) 

where  F( • ) is  the  Fourier  transform  of  the  processor  output 
density  function.  The  direct  computation  of  the  transform 
of  the  expression  given  by  Eq  (32)  is  a mathematical  exer- 
cise beyond  the  scope  of  this  thesis.  A search  of  the 
Campbell  and  Foster  transform  table  (Ref  1)  failed  to  iden- 
tify the  necessary  transform  pair,  which  may  indicate  that 
a general  transform  for  the  probability  density  function 
given  in  Eq  (32)  does  not  exist.  The  lack  of  a suitable 
transform  with  which  to  obtain  the  non-linear  approximate 
processor  output  characteristic  function  means  that  analytic 
performance  expressions  for  the  general  non-linear  approx- 
imate processor  cannot  be  derived.  Single  channel  perform- 
ance expressions  can  be  defined  by 

00 

w 

I’fa=  £ <35) 

but  the  required  integration  is  very  complex  and  its  computa- 
tion viewed  as  of  little  real  value. 

If  it  is  assumed  that  the  conditional  detector  output 
variances  are  equal,  then  b^  = 0 and  the  single  channel 
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processor  output  density  function  is  then  given  by 


fz(z) 


|Z| 

(ci-ei)2 


exp 


X I, 


(ci-ei) 


(36) 


While  this  expression  is  less  complex  than  that  given  by 
Eq  (32) » a general  Fourier  transform  pair  still  could  not  be 
found  and  the  direct  integration  of  the  expression  remains 
overly  complex.  Here  again  the  complexity  of  the  necessary 
mathematics  prevents  the  acquisition  of  analytic  performance 
expressions  for  the  non-linear  approximate  processor.  The 
attempt  to  find  analytic  performance  expressions  is  next 
directed  at  the  simpler  case  of  the  linear  approximate 
processor. 

Linear  Approximate  Processor 

Beginning  with  the  linear  approximate  processor  algorithm 
given  by  Eq  (24),  the  processor  output  function  for  a single 
channel  "i"  is  given  by 

z(y)  = [bi  + i(ci2-ei2)]  y (37) 

Through  transformation  of  variables,  the  unconditional  single 
channel  processor  output  probability  density  function  is 
expressed  by 
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the  range  of  z depending  upon  Ki,  where  Ki  = + ^(c^^-e^^). 

Using  transform  pair  655*1  from  the  Campbell  and  Foster 
transform  table  (Ref  1),  the  Fourier  transform  for  an  M 
channel,  N look  linear  approximate  processor  is  given  by 
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where  j is  the  imaginary  unit  . Unfortunately,  an  inverse 
transform  for  Eq  (39),  which  would  yield  an  expression  for 
the  general  multichannel  linear  approximate  processor  out- 
put joint  density  function,  could  not  be  found  except  for 
the  special  case  where  each  channel  is  assumed  to  be  iden- 
tically distributed.  Under  this  assumption  Eq  (39)  can  be 
reduced  to 


F(f) 


exp 


N M m 

- N M yj 

20^ 

exp 

2a^(  1+ j4TTKa^f ) 

(l  + j4TTKa^f) 


(40) 


where  the  density  parameters  are  no  longer  channel  dependent. 
Using  transform  pair  650.O  from  the  Campbell  and  Foster  table, 
the  processor  output  joint  density  function  for  M identically 
distributed  channels  is  given  by 
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While  the  joint  density  function  is  in  analytic  form,  the 
calculation  of  the  performsince  parameters  Pq  and  requires 
that  this  pdf  be  integrated  according  to  Eqs  (3^)  (35)  • 

This  integration  is  prohibitively  complex  and  a search  of 
available  integral  tables  failed  to  yield  an  integral  of  the 
required  form.  Even  for  the  simpler  structure  of  the  linear 
approximate  processor,  the  derivation  of  exact  performance 
expressions  remains  prohibitively  cumbersome.  When  exact 
performance  calculations  are  impossible,  it  is  often  useful 
to  determine  bounds  on  error  probabilities. 

Chemoff  Bound 

An  exponetially  tight  bovind  on  the  probability  of  false 
alarm  and  the  probability  of  detection  is  the  Chernoff  bound 
(Refs  2:126;  12:121).  This  bound  is  obtained  from  the  moment 
generating  function,  }a(s),  which  is  derived  from  the  condi- 
tional characteristic  function  ^(p|Ho).  The  moment  generating 
function  and  the  resultaint  bounds  on  the  probability  of  false 
alarm  and  the  probability  of  detection  are  given  by 

}i(s)  = ln[E(e®^|Ho)]  = In  $(-js|Ho)  (42) 
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Pj3  - 1-exp  [^(s)  + (l-s)fi(s)] 
PpA-  exp  [}i(  s ) - S|a(  s )] 


(43) 


(44) 


where  E(*  ) indicates  expected  value  and  p(s)  = ^|i(s)/<)s  . 
Using  Eqs  (33)  and  (39) i the  characteristic  function  for  the 
linear  approximate  processor  can  he  derived.  Because  no 
transform  for  the  non-linear  approximate  processor  output 
density  function  exists,  neither  the  moment  generating  func- 
tion or  the  Chernoff  bound  for  the  non-linear  approximate 
processor  can  be  obtained. 

The  conditional  characteristic  function  for  the  linear 
approximate  processor  output  is  given  by 


KhIHq)  = 


/ ®i0 

exp ^ I expl 

M \ ZUig 

n — ^ — 

i=l 


m 


iO 

2ai§(l-j2KiGi§p) 


(l-j2Kiai0H) 


N 


(45) 


for  ]i  real,  from  which  the  moment  generating  function 


M 

u(s)  = N S 
i=l 


niiO 


L 2Ui§  \ l-2Kiai§s 


- 1 I - ln(l-2Kiai(^s) 


(46) 


is  obtained.  The  partial  derivative,  ]a(s),  is  found  to  be 


M 

u(s)  = N Z 
i=l 


2KiCTi§  “iO^i 


1_  l-2KiaioS  (l-2Kiaios)^. 


(4?) 
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The  bounds  on  Pq  and  Pp^  are  now  defined  by  Eqs  (43)  aind  (44), 
where  the  optimal  value  of  s is  the  solution  to  the  equation 

lj(s)  = w 

where  w is  the  previously  defined  threshold. 

The  use  of  the  bounds  in  general  require  that  a poly- 
nomial be  solved  to  obtain  the  optimal  value  of  s.  In  all 
but  the  most  simple  cases,  this  can  be  quite  cumbersome. 

For  the  single  channel,  one  look  case,  the  optimal  value  of 
s was  determined  to  be 


s 


(w-KiQi^)  + ^Ki^aio+mioKiW ' 
2Kiai§w 


(49) 


Assuming  for  simplicity  that  m^Q  = 0 eind  that  the  variances 

p 

under  each  hypothesis  are  equal,  s = 2/rij^  where 

and  the  bound  on  the  probability  of  false  alarm  is  given  by 


^FA 


21n(v)+ni 


(50) 


where  v is  the  linear  threshold.  The  limit  of  this  bound  as 
r|j^  increases  without  bound  is 

PpA-  1/e  = 0.37  (51) 

Intuitively,  one  would  expect  the  bound  on  Pp^  to  approach 
zero  as  increased  without  bound,  since  the  conditional 
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processor  output  density  functions  separate  as  increases. 
The  loose  bound  on  Pp^  of  O.37  indicates  that  the  Chernoff 
bound  obtained  is  valid  only  for  small  values  of  , a result 
consistent  with  the  small  argument  assumption  originally  used 
in  the  derivation  of  the  linear  approximate  processor. 

Numerical  Analysis  , 

The  lack  of  analytic  expressions  with  which  to  compare 
between  signal  processors  necessitated  the  use  of  numerical 
methods  to  achieve  a common  basis  for  comparison.  The  pro- 
cedure used  here  was  to  write  a FORTRAN  computer  prograjn 
which  directly  computed  the  conditional  processor  output 
densities  of  each  processor  for  a given  set  of  parameters. 
These  densities  were  then  numerically  integrated  over  a 
range  of  threshold  values  to  generate  an  array  of  values  for 
the  probability  of  detection  and  the  probability  of  false 
alarm.  These  values  were  in  turn  plotted  as  a receiver 
operating  curve  ( ROC ) . 

The  computation  of  joint  density  functions  for  multi- 
channel receiver  performance  evaluation  made  use  of  a fast 
Fourier  transform  (FFT)  algorithm  to  compute  the  transforms 
of  the  conditional  processor  output  densities  for  each  chan- 
nel. The  product  of  the  individual  channel  transforms  was 
then  inverse  transformed  to  obtain  the  conditional  joint 
processor  output  density  functions.  These  joint  density 
functions  were  then  numerically  integrated  as  above  to  obtain 
values  for  Pjj  and  Pp^  which  were  plotted  as  a receiver  op- 
erating curve . 


The  conditional  output  density  functions  for  the  non- 
linear approximate  processor  and  the  linear  approximate  pro- 
cessor were  derived  from  Eqs  (32) i (36) » and  (38)  and  were 
expressed  in  terms  of  the  parameters  = mj^o/a^^Q,  = 
mii/Qii , and  p = . Single  channel  processor  perform- 

ance was  calculated  for  changes  in  target  hypothesis  mean 
and  for  changes  in  target  hypothesis  variance.  The  param- 
eters  used  in  this  analysis  were  chosen  primarily  for  con- 
vienience  in  computation  and  are  scaled  versions  of  the 
pareimeters  estimated  from  the  Lockheed  Background  Measure- 
ments Program  data  discussed  in  Chapter  II.  Since  the  values 
2 

of  m^Q  and  estimated  from  the  Lockheed  data  are  linearly 
related  to  the  detector  output  voltage  (Ref  6:3-27),  it  was 
asserted  that  the  ratio  of  detector  output  mean  to  variance 
(tiq)  could  effectively  be  set  for  a given  value  through 
proper  selection  of  processor  components.  The  parameters 
used  to  observe  the  effects  of  changes  in  mean  upon  processor 
performance  are  listed  in  Table  II.  The  parameters  used  to 
observe  the  effects  of  changes  in  variance  upon  processor 
performance  are  listed  in  Table  III. 


Table  II 

Channel  Parsimeters  for  Change  in  Mean 


^0 

P 

b.O 

T^To 

1 .0 

8.0 

24.0 

1 .0 

8.0 

32.0 

1.0 
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Table  III 

Channel  Parameters  for  Change  in  Variance 


^0 

m 

p 

8.0 

WTo 

1 .0 

8.0 

24.0 

1.5 

8.0 

24.0 

2.0 

Typical  non-linear  processor  output  densities  for  change 
in  mean  are  illustrated  in  Figures  12-  14.  It  is  easily  ob- 
served that  the  conditional  processor  output  density  func- 
tions separate  as  the  target  hypothesis  meein  increases.  The 
resulting  receiver  performance  curves  are  shown  in  Figure  15. 
For  a probability  of  false  alarm  equal  to  0.10,  the  proba- 
bility of  detection  increases  from  0.45  to  O.96  as  ri^  goes 
from  16.0  to  32.0  with  r]Q  equal  to  8.0  . Typical  non-linear 
processor  output  densities  for  change  in  target  hypothesis 
variance  are  illustrated  in  Figures  I6  and  17.  The  resulting 
pe". ormance  curves  are  shown  in  Figure  18.  As  the  target 
hypothesis  variance  increases  to  twice  the  null  hypothesis 
variance  with  constant  mean  to  variance  ratios,  the  proba- 
bility of  detection  goes  from  0.80  to  0.99  for  a probability 
of  false  alarm  of  0.10  . The  observed  changes  in  performaince 
indicate  that  the  receiver  is  more  sensitive  to  changes  in 
mean  than  to  changes  in  variance.  A representative  FORTRAN 
program  used  to  generate  the  ROC ' s for  the  non-linear  approx- 
imate processor  is  listed  in  Appendix  B. 

To  determine  the  effect  of  multiple  channels  upon  re- 
ceiver performance,  performance  curves  were  generated  for  one, 
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Fig.l2.  Non-linear  Processor  Output  Densities:  r|Q=8.0, 
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Fig. 13*  Non-linear  Processor  Output  Densities:  rio=8.0,  ni  = 24. 


Fig.  l4.  Non-linear  Processor  Output  Densities:  no=8*0»  Di=32. 


Fig.  16.  Non-linear  Processor  Output  Densities:  no=8.0,  nj^  = 24. 
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Fig,  17.  Non-linear  Processor  Output  Densitiesj  ri„  = 8.0,  ri.=24. 


two,  ■ nd  three  channel  receivers.  Three  nearly  identical 
channels  were  used;  the  channel  parameters  are  listed  in 
Table  IV.  The  performance  curves  for  the  individual  channels 
are  plotted  in  Figure  19*  For  a probability  of  false  alarm 
of  0.10,  each  of  the  individual  channels  has  a probability 
of  detection  of  approximately  0,49.  Multichannel  receiver 
performance  curves  for  the  non-linear  approximate  processor 
are  shown  in  Figure  20.  The  probability  of  detection  in- 
creases to  0.66  for  two  channels  and  to  0.82  for  three  chan- 
nels, ax\  improvement  of  67%  over  single  channel  performance. 
Since  a similar  joint  density  function  would  be  obtained  for 
a multi-look  receiver,  by  extension  of  the  multichannel 
results,  one  would  also  expect  receiver  performance  to  im- 
prove for  a multi-look  receiver.  A representative  FORTRAN 
program  used  to  generate  multichannel  ROC ' s is  listed  in 
Appendix  C . 

The  saime  procedures  were  applied  to  the  linear  approxi- 
mate processor  using  identical  parameters.  Typical  linear 
processor  output  densities  for  change  in  mean  are  illustrated 
in  Figures  21  -23.  The  resulting  receiver  performance  curves 
are  shown  in  Figure  24,  The  observed  performance  is  identi- 


Table  IV 

Multichannel  Parameters 


Channel 

’^O 

’^l 

P 

1 

8.0 

16.0 

1.0 

2 

9.0 

18.0 

1 .0 

1 

10.0 

20.0 

1 .0 
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Fig.  20.  Multichannel  Non-linear  Processor  ROCs 
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Fig.  22.  Linear  Processor  Output  Densitiesi  r|Q=8.0,  r\^  = Zk.O,  p=l . 
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Fig.  23.  Linear  Processor  Output  Densities;  r|„  = 8.0,  r).=32. 


cal  to  that  of  the  non-linear  approximate  processor.  Typical 
linear  processor  output  densities  for  change  in  variance  are 
illustrated  in  Figures  25  and  26.  The  resulting  performance 
curves  are  shown  in  Figure  2?.  Again,  the  performance  of 
the  linear  approximate  processor  is  identical  to  the  perform- 
ance of  the  non-linear  approximate  processor.  A representa- 
tive FORTRAN  program  used  to  generate  the  ROC ' s for  the 
linear  approximate  processor  is  listed  in  Appendix  D.  Multi- 
channel receiver  performance  durves  for  the  linear  approximate 
processor  are  shown  in  Figure  28.  Linear  and  non-linear  pro- 
cessor multichannel  performance  is  also  identical.  The  equal 
performance  result  is  proven  analytically  for  the  single 
channel  case  in  Appendix  E. 

Ad-hoc  Linear  Processors 

An  ad-hoc  signal  processor  is  a processor  whose  struc- 
ture is  based  upon  engineering  intuition  rather  than  statis- 
tical analysis.  The  ad-hoc  linear  processors  presented  in 
Chapter  III  are  used  here  to  compare  the  relatively  simple 
ad-hoc  processor  structure  with  the  performance  of  the  non- 
linear and  linear  approximate  processors  derived  from  signal 
statistics . 

The  conditional  processor  output  density  functions  for 
the  ad-hoc  linear  processor  #1  are  obtained  from  the  processor 
algorithm  given  by  Eq  (25)  through  transformation  of  vari- 
ables. The  resulting  single  channel  conditional  processor 
output  density  fimctions  are  given  by 
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Fig.  25.  Linear  Processor  Output  Densities;  ^1=24.0, 


RECEIVER  OPERATING  CURVE 


o 

o 


PROBABILITY  OF  FALSE  ALARM 


Fig.  28.  Multichannel  Linear  Processor  ROCs 
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where 

These  density  functions  can  be  expressed  in  terras  of  the 
parameters  iIq,  , and  p previously  defined.  Since  the  pro- 
cessor structure  is  similar  to  that  of  the  linear  approxima~ 
processor,  differing  only  in  the  value  of  K^,  performance 
curves  for  the  ad-hoc  linear  processor  were  generated  by 
using  the  program  listed  in  Appendix  D with  only  slight  mod- 
ification. The  resulting  performance  curves  for  the  ad-hoc 
linear  processor  #1  for  changes  in  mean  and  variance  are 
shown  in  Figures  29  and  30-  For  identical  parameters,  the 
performance  of  the  ad-hoc  linear  processor  #1  is  identical  to 
the  performance  observed  for  both  the  non-linear  and  linear 
approximate  processors. 

Similarly,  the  single  channel  conditional  processor  out- 
put density  functions  for  the  ad-hoc  linear  processor  #2  were 
derived  from  the  processor  algorithm  given  by  Eq  (26).  The 
density  functions  for  the  ad-hoc  linear  processor  #2  are  of 
the  same  form  as  Eqs  (52)  and  (53)  except  that  Kj^  is  now 
defined  by 
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These  density  functions  can  he  expressed  in  terms  of  the 
parameters  ni^Q,  , and  p.  For  the  purpose  of  performance 
calculations  the  value  of  m^^Q  was  arbitrarily  chosen  to  be 
2.5  . To  make  the  performance  calculations  for  the  ad-hoc 
linear  processor  #2  consistent  with  those  for  the  other 
processors,  m^^j^  was  determined  by  the  following  equation* 

""il  = "“ioP  ^ ^54) 

where  r^Q,  , and  p are  the  channel  parameters  listed  in 
Tables  II  and  III.  As  before,  the  program  listed  in  Appendix 
D was  modified  to  generate  the  ROC's  for  this  processor.  The 
resulting  performance  curves  for  the  ad-hoc  linear  processor 
#2  for  changes  in  mean  and  variance  are  shown  in  Figures  31 
and  32.  For  identical  parameters,  the  performance  of  this 
processor  is  also  identical  to  the  performance  observed  for 
both  the  non-linear  and  linear  approximate  processors. 

Performance  Parameter  Dependence 

In  an  effort  to  determine  the  manner  in  which  receiver 
performance  depends  upon  the  parameters  no»  » and  p, 
functions  of  these  parameters  were  tested  against  previous 
results. 

The  first  function  tested  was  parameters 

listed  in  Table  II,  this  function  has  the  values  8.0,  I6.0, 
and  24.0  . To  test  the  dependence  of  receiver  performance 
upon  this  function,  was  set  equal  to  0.0  and  set  equal 
to  8.0,  16.0,  and  24.0  . The  resulting  performance  curves 
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for  the  non-linear  approximate  processor  are  shovm  in  Figure 
33.  Comparing  performance  with  that  of  Figure  I5 , it  is 
obvious  that  receiver  performance  does  not  depend  upon  the 
difference  between  the  hypothesis  meain  to  variance  ratios. 

The  next  function  tested  after  analyzing  the  previous 
results  was  • For  the  parameters  listed  in  Table 

II,  this  function  has  the  values  1.17,  2.07,  and  2.82  . To 
test  the  dependence  of  receiver  performaince  upon  this  func- 
tion, riQ  was  set  equal  to  6.0  and  set  to  13.112,  20.^3, 
and  27.856  . The  resulting  performance  curves  for  the  non- 
linear approximate  processor  are  shown  in  Figure  3^.  Compar- 
ing these  performance  curves  with  those  of  Figure  I5 , it  is 
apparent  that  performance  is  identical.  One  may  conclude 
from  this  result  that  receiver  performance  depends  upon  the 
difference  between  the  square  roots  of  the  hypothesis  mean 
to  variance  ratios.  This  differs  from  the  Gaussian  signal 
case  where  receiver  performance  depends  upon  the  difference 
between  the  hypothesis  mean  to  standard  deviation  ratios, 
^il^'^il  “ "^iO'^'^iO  • Determination  of  receiver  performance 
dependence  upon  functions  also  involving  p were  unsuccessful. 
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Fig*  33*  Functional  Receiver  Performance  Dependence, 
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Conclusions  and  Recommendations 


The  results  of  this  thesis  indicate  a number  of  conclu- 
sions and  also  suggest  a number  of  areas  for  further  study. 

Conclusions 

The  following  conclusions  are  made  by  this  study: 

The  non-central  chi-square  statistical  detector  output 
model  appears  to  be  valid.  The  model  has  a realistic  physical 
basis  in  the  complex  Gaussian  random  process  used  to  model  the 
received  infrared  field  and  is  further  supported  by  available 
experimental  infrared  background  statistical  data. 

Receiver  performance  equal  to  that  of  either  the  non- 
linear approximate  processor  or  the  linear  approximate  pro- 
cessor derived  from  the  optimal  signal  processor  can  be 
obtained  using  a relatively  simple  ad-hoc  linear  reciever. 

Receiver  performance  depends  upon  the  difference  between 
the  square  roots  of  the  hypothesis  mean  to  variance  ratios 
(V^il/Gii  W^iO^’^iO^  ratio  of  hypothesis  variances 

(Oif/crJ). 

Receiver  performance  improves  as  the  number  of  receiver 
channels  and  the  number  of  looks  increase . 


Recommendations 

The  following  are  recommendations  of  areas  for  further 
study  pertaining  to  multichannel  infrared  receiver  perform- 
ance: 

While  the  statistical  detector  output  model  appears  to 
be  valid,  raw  experimental  infrared  background  data  should  be 
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ainalyzed  to  validate  the  model  for  different  types  of  back- 
grounds. This  might  be  accomplished  by  analyzing  the  raw 
LMSC  Background  Measurements  Program  data  or  by  original 
experiment . 

Another  area  for  futher  study  is  to  determine  a statis- 
tical model  for  the  infrared  emissions  of  anti-aircraft 
missiles  and  to  determine  the  effects  of  this  model  upon  the 
overall  signal  statistics. 

The  attempts  to  obtain  analytic  performance  expressions 
were  frustrated  by  the  complexity  of  the  mathematics  involved 
with  the  non-central  chi-square  probability  density  function. 
It  may  be  possible  to  simplify  the  required  mathematics  and 
gain  some  insight  into  receiver  performance  dependence  by 
using  Gaussian  approximations  to  the  non-central  chi-square 
pdf.  Possible  approximations  to  the  non-central  chi-square 
pdf  are  suggested  in  a University  of  Michigan  technical  re- 
port titled  "Approximations  to  the  Noncentral  Chi  Square 
Distributions  with  Applications  to  Signal  Detection  Models" 
(Ref  9) . 

While  it  has  been  determined  that  the  four  signal  pro- 
cessor structures  examined  in  this  thesis  have  equal  perform- 
ance, it  is  not  known  how  this  performance  compares  with  that 
of  the  optimal  processor.  It  is  possible  to  calculate  opti- 
mal processor  performance  numerically  by  first  computing  the 
conditional  detector  output  density  functions  and  then  process 
these  density  functions  according  to  the  optimal  processor 
algorithm  to  obtain  the  conditional  processor  output  density 
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functions.  These  pdf's  can  then  be  numerically  integrated 
to  obtain  receiver  operating  curves  comparable  to  those  of 
the  other  processors. 

Other  areas  for  futher  study  might  include  temporal 
signal  processing,  adaptive  receivers,  and  determination  of 
parameter  dependence  for  multichannel  and  multi-look  receiver 
performance . 
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Appendix  A 

The  combined  relative  frequency  plots  for  filters  4,  5i 
and  6 are  illustrated  in  Figures  A-1 , A-2,  and  A-3*  These 
plots  closely  approximate  the  experimental  histograms  obtained 
by  LMSC  for  the  corresponding  filter.  Figures  A-4  and  A-5 
are  plots  of  the  estimated  pdf  parameters  as  a function  of 
filter  bandwidth.  There  appears  to  be  a nearly  lineair 
relationship  between  density  parameters  and  filter  bandwidth, 
with  deviations  resulting  from  parameter  estimation  error. 

For  Figures  A-4  and  A-5,  the  left  edge  of  each  filter  is 
fixed  at  4.39  V™  in  all  cases. 
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Fig.  A-4.  Background  Mean  Values 
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Fig.  A-5.  Background  Standard  Deviation 
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Appendix  B 


The  following  program  listing  is  typical  of  the  FORTRAN 
programs  used  to  generate  receiver  performance  curves  for  the 
non-linear  approximate  processor. 


PROGRAM  GRAPH  ( IN  PUT , OUTPUT , PLOT) 

IHTFGFR  POINT 

DIMENSION  Y(602)»03(502), 01(60?) 

DIMENSION  PFA (40) ,oD (40) 

DATA  FAMIN,DMIN,-An, DO/0 . 0 , 0 . 0,  0 . 25,  0 . 16657/ 

CALL  PLOT(0.,-3.0 ,-3) 

CALL  PLOT(0.,2.5, -3) 

CALL  PLOT(0. 0,6.0  ,3) 

CALL  PLOT(4.0,6.0,2) 

CALL  PLOT(4.0,0.0 ,2) 

CALL  PLOT(4.0,6.0 ,3) 

CALL  PLOKO.O,  0.0 ,2) 

CALL  AXIS(O.,O.,26HPPORA0ILITY  OC’  FALSE  ALARM, -26 , 4 . 0 , 

♦ 0.0,«=‘AMIN,FAO) 

CALL  AXIS(0.,0.,24^PROPARILITY  0-  DETECT I3N , 24, 6 . 0 , 90 . , 
♦OMIN,nn) 

CALL  SYMBOL(0.25,6.4,0.14,24MRECEIVER  OPERATING  CURVE, 
4-0.0, 24) 

* POINT  = NUMBER  OF  DATA  POINTS  IN  EACH  DENSITY 

* INDEX  J = NUMBER  0^  DATA  SETS 

* 

POINT=240 

ETA0=8.0 

ETA1=24.0 

RHO=1.0 

* 

* COMPUTE  HO  AND  HI  DENSITIES  FOR  EQUAL  VARIANCE 

OIO=?.0*((ETA1/RHO)-2.0*SORT (ETA1*ETAO/RHO) ♦ETAO) 

CIO  = SQRT(ETAO)/(SQRT  (ETAl/RHO)  -SDRKETAO)  ) 

DI1=RHO»OIO 

CI1=S0RT(ETA1)/ (SORT (ETAl) -SORT ( ETAO* RHO)  ) 
IF(CIO.LT.O.O)  CI0=ABS(CI0) 

IF(CIl.LT.O.O)  CI1=ABS(CI1) 

DO  50  1=1, POINT 
YCI)=0. l*(FLOAT(I-l)  ) 

DO (I)=XCHINL (Y(I)  , ETAO, CIO, DIO) 

01(I)=XCHINL(Y(I)  ,ETA1,CI1,0I1) 

50  CONTINUE 

* 

* calculate  roc  points 

* 

N=0 

DO  75  M = 30,138,4 
N=N41 

PFA (N)=GRAL (DO ,M, POINT) 

PO(N) =GRAL{ 01, M, POINT) 

75  CONTINUE 
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* 

* 


PLOT  THE  RECEIVED  OPERATION  OIJPV-  (ROC)  » PO  VS  PFA 


PFA (N»l) =PO(N+l) =0.0 

PFA()U2)=0.25 

«»0(N*2)  =0.16667 

CALL  LINE(PFA,PO,N, 1,5,1) 

* 

♦ CONTINUE  FOR  DENSITIES  WITm  UNEQUAL  VARIANCE 

* 

read*, POINT 
DO  300  J=2,3 
REAn‘,ETAO,ETAl,RHO 

* 

* COMPUTE  HO  AND  HI  DENSITIE? 

* 

A=SnPT{ETAl) -SORT (ETAO^RHO) 

B=SQRT(ETA1/RH0) -SORT (FTAO) 

C=l.C-( l.O/RHO) 

O=RHO-1.0 

E=(E'!’A1/RH0)  ♦ETA3-(2.0*SnR'r(ETAQ»ETAl/RH0)  ) 

^=ETAH-(RHO*ETAO)  -(2.0*SQRT(ETAO»ETA1»RHO)  ) 

3=2.0»E/C 

00  ICO  1=1, POINT 

V{I)=0. 1*FL0AT (I-l) 

!)0(I)=XNL(Y(I)  ,ET  AO  , RHO, A ,C  ,F  ,r,  ,9  ,E) 

Dl(I)  =XNL  (Yd)  ,Er  A1,RH0,A,0,'^,S,  A,F) 

100  CONTINUE 

CALL  SCALE  (00,4.3 , POINT, 1) 

CALL  SCALE(01,4.0 ,POINT,l) 

THR=O0(POINTi-2) 

« 

• FIND  THE  STARTING  OOINT  OF  HO  DENSITY , INDEX  K 

* 

DO  110  1=1, POINT 
IF(O0 (I) .GT.THR)  GO  TO  120 
110  CONTINUE 

120  K=(INT(FL0AT (I)/2.0) )*2 

• FIND  END  POINT  0-  HI  DENSITY,  INDEX  L 

« 

THR=ni (POINT+2) 
no  130  1=1, POINT 
IF(D1(I)  .GT.THR)  GO  TO  140 
130  CONTINUE 
140  DO  150  L=I, POINT 

IF(ni (L)  .LT.THR)  GO  TO  160 
150  CONTINUE 

160  L=(INT(FLOAT(L)/2.0) ) *2 
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THIS  SFCTION  DIVIDES  THE  nT'^PE'»ENCE  BETWEEN  THE  CHOSEN 
BEGINNING  AND  END  ^^OINTS  OE  THE  DENSITIES  BY  30  TO  GET 
THE  LOOP  INCREMENT  INDEX,  I,  SO  AS  TO  CALCULATE  30  ROC 
POIN'^S  OVER  THE  SIGNIFICANT  PORTIONS  OF  THE  DENSITIES 

I*IN^ (FLOAT (L-K) /30 . ♦0.5) 

IF(I.EO.O)  1=1 

PRINTS, ’'ROC  PARAMETERS  **,K,L,T 
PRINT*,  ••  •• 

•4=0 

BO  200  M=<,L,I 
N=N^1 

PFA(N)=GRAL(D0,M, POINT) 

Pn(N)=GRAL (01,M,30INT) 

200  CONTINUE 

PFA (N^1)=PO(N+1)=0.0 
PFA(N*2)=0.25 
PD(N+2) =0.15667 

* 

♦ PLOT  THE  RECEIVER  OPERATION  C'JPVF  (ROC)  I PO  VS  PFA 

♦ 

CALL  LINE(PFA,P0,N,1,5, J) 

300  CONTINUE 

CALL  PLOTE(N) 

STOP 

END 


THIS  FUNCTION  COMPUTES  THE  DUTpJT  DENSITY  EOR  THE  NON- 
LINEAR PROCESSOR  WHEN  THE  VAOIANCES  ARE  NOT  EQUAL  FOR 
THE  TWO  HYPOTHESES,  IE.  RHO  IS  GREATER  THAN  ONE. 

FUNCTION  XNL( Y,ETA,R, AA,B3,CC,nD,EE,FF) 

DOUBLE  V 

V=(Y4 (2.»00)-2.*SDRT (Y*OD^DD*DD) ) /BB 
V=-(V+(ETA/2.0) ) 

X=(1./BB)-(AA/(BB*SDRT{(Y*(R-1.)  /2.0)  ♦CO)  ) 

X=ARS(X) 

Z1=X*DEXP(V) 

Z2=2.0*SQRT(ETA)rBB 

Z2=Z?* (SQRT( ( Y*BB/2. 0) ♦FF) -pc) 

Z2=FI0 (Z2) 

XNL=Z1*Z2 

RETURN 

END 
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♦ THIS  FUNCTION  CO'^PUTES  A POINT  ON  THE  NON-LINEAR 

♦ PROCESSOR  OUTPUT  DENSITY  i<HEN  VARIANCES  ARE  EQUAL 

♦ 

FUNCTION  XCHINL(Y  ,ETA,CI,OI) 

OOUBLF  V 

V=-((Y**2) /OI+ETA/2. 0) 

Zls2.0*OEXP(V)*A3S(Y)/DI 

72=FI0{Y»CI) 

XCHINL=Z1*Z2 

RETU'^N 

END 


* THIS  FUNCTION  COMPUTES  THE  MODIFIED  BESSEL  FUNCTION  OF 

* THE  FIRST  KIND  O*^  ORDER  ZERO 

* 

FUNCTION  FIO(X) 

DIMENSION  XL(6),XH(8) 

DOUBLE  W 

DATA  XL/3.5156229 ,3. 08DP424,1.2057492,0.2S59732, 0 .0360763, 
10.  cot* 13/,  XH/ 0.0  132  8 592,  0.  0 0 2?5319,- 0.0  01 57565,  0.  0 0916281, 
2-0.02057706,0.02635537,-0.01 6476^3,0. 00392377/ 

H=X 

XA=ABS(X) 

IF  (XA.GT. 0.0001)  SO  TO  10 

FI0=1.0 

RETURN 

10  IF  (X.GT.O)  GO  TD  20 
X=XA 

20  T=X/3.75 

IF  (T.GT.1.0)  GO  TO  40 

FI0=1.0 

DO  30  J=l,6 

FI0=‘^I0+XL(J)*(T**(2*J)) 

30  CONTINUE 
RETUON 

40  FI0=0. 39894228 
TI=1.0/T 
DO  50  J=l,8 
FI0=FT0+XH(J)*(TI»*J) 

50  CONTINUE 

FI0  = FI0 *OEXP(W)* (1.0 /SORT (X)  ) 

RETURN 

END 
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THIS  FUNCTION  COMPUTES  THE  iNTECRfiL  OF  THE  FUNCTION 
CONTaiNEO  IN  THE  ARRAY  A USINS  A SIMPLE  RECTANGULAR 
SUMMATION  ALGORITHM 

FUNCTION  GRAL  (A,J,K) 

DIMENSION  A(K) 

M = K-? 

GRAL=0.0 
00  100  I=J,M,2 
GRAL=GRAL»A(I)+A( 1+2) 

100  CONTINUE 

GRAL=0. l^GRAL 
RETURN 

end 
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Appendix  C 

The  following  program  listing  is  typical  of  the  FORTRAN 
programs  used  to  generate  multichannel  receiver  performance 
curves  for  the  non-linear  and  linear  approximate  processors. 
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PROGRAM  GRAPH  (I NPUT, OUTPUT, PLOT) 

INTEGER  POINT  _ _ 

COMPLEX  A ( 512) ,8(512) ,C( 512) ,0(512) 

DIMENSION  r(51U) , DO (514) ,01(514)  , 

DIMENSION  PFA(50 ) ,PO(50) 

DATA  FAMIN, DMIN, FAD, DD/0. 0,0  .0,0. 25,0. 16667/^  ... 

CALL  PLOT(0..-3.0,-3) 

CALL  PLOT(0.,2.5,-3)  

CALL  PLOT(O.Q,6.a,3) 

CALL  PLOT( 4.0,6.0 ,2)  _ _ 

CALL  PLOT(4.0, 0.0,2) 

CALL  PLOT(4.0,6.Ct3) 

CALL  PLOT(0.0,0.0,2) 

CALL  AXIS( 0 ,,0.,26HPRO8A8ILITY  OF  FALSE  AL ARM,-2 6,4. 0 t 
♦ 0.0.FAMIN,FAD) 

CALL  AXIS(0.,0.,24HPROBA3ILITY  OF  DETECTION, 24, 6 .0 ,90 . i 
'«>0HIN,00) 

CALL  SYMBOL(0.25,6.4,0.14,24HRECEIVER  OPERATING  CURVE, 
4-0.0,24) 

POINT  = NUMBER  OF  DATA  POINTS  IN  EACH  DENSITY 

INDEX  J = NUMBER  OF  DATA  SETS  _ 

IFFT  = NUMBER  OF  FFT  POINTS 

NFFT  s LOG2(IFFT)  

NCHAN  = NUMBER  OF  CHANNELS 

DELTA  = DIFFERENCE  BETWEEN  DENSITY  POINTS  

DO  300  J=l,3  . . „ 

READ*, POINT , IFFT , NFFT, NCHAN 

INITIALIZE  COMPLEX  MULTIPLIER 

00  10  1=1, IFFT 

B(I)  = (l,  0,0.0)  _ „ 

0(I)  = (1.0,0  .0) 

10  CONTINUE  _ _ _ _ 

OELTA=0.1 

00  50  JJ=1, NCHAN  __  __ 

READ*,ETAO,ETA1,RHO 

COMPUTE  CHANNEL  DENSITIES  

DIC=2,0* ((ETAl/RHO) -2.0* SORT (ETA1*ETAO/RHO) fETAO) 

CIC=SQRT(ETAO)/(SQRT(ETA1/RHO) -SQRT(ETAO))  _ 

DI1=PHO*DIO 

CI1=S0PT(ET Al) /(SQRT(ETAl) -SORT ( ETAO* PHO) ) 

DO  20  1=1, POINT 

YII)=OELTA*FLOAT(I-l)  _ 

ACI)=XCHINL (Y(I) ,ETAO,CIO,OIO) 

C<  I)=XCHINL(Y(I)  ,ETA1,CI1,0I1) 

20  CONTINUE 
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♦ ZERO  FILL 

« 

KK=POINT^l  ■ 

00  30  I=KK,IFFT 
Y(I)=DELTA*FL0AT(I-1» 

AfI)-(O.OtO  .0) 

C(I)  = (0.0«a.0) 

30  CONTINUE 

♦ 

* TAKE  FFT  OF  DENSITIES 

« 

CALL  FFT(A,NFFT,IFFT,-1) 

CALL  FFr(C,NFFT,IFFT,-l) 

« 

* PERFORM  CONVOLUTION  BY  MULTIPLYING  TRANSFORMS 

DO  40  I=1,IFFT 
a<i»=A(i»*9(i) 

D(I>=C(I)*D(I>  . 

40  CONTINUE 

* 

• RETURN  FOR  NEXT  CHANNEL 

50  CONTINUE  ' ' ’ 


• TAKE  INVERSE  TRANSFORM  TO  GET  JOINT  DENSITIES 


CALL  FFT{B,NFFT»  IFFT,«-1» 

CALL  FFT(0,NFFT,IFFT,*1)  . ^ 

* SCALE  AND  SET  TO  REAL  VALUES 

• T SCALES  DATA  TO  ACCOUNT  FOR  DIFFERENCES  IN  THE  DELTAS 

♦ SINCE  THE  FFT  ALGORITHM  ASSUMES  A CELTA  OF  1 

T*OELTA**NCHAN 

00  60  I=ltIFFT  _ _ __  _ 

OO  (I)=T*B<  D/IFFT 

01(I)=T*D(  D/IFFT  . _ 

60  CONTINUE 

♦ COMPUTE  ROC  FOR  MULTICHANNEL  PROCESSOR 

• DENSITY  IS  NOW  COMPOSED  OF  IFFT  POINTS  

» 

POINT=IFFT  . _ 

CALL  SCALE(O0,4.0»POINT»l) 

CALL  SC4LE(Dl,4.0fPOINTfl) 

THR=0CCP0INTf2) 
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* 

« 

» 


FIND  THE  STARTING  POINT  OF  HO  DENSITY , INDEX  K 
DO  lie  1=1, POINT 

IF(Dfl(I)  .GT.THR)  GO  TO  120  

110  CONTINUE 

120  K=(INT(FLOAT (I)/2.0) 

IF(K.LE.O)  K=1 

FIND  END  POINT  OF  HI  DENSITY,  INDEX  L 
THR=Dl(POINT*2) 

DO  13C  1 = 1, POINT  _ . 

IF(01(I»  .GT.THR)  GO  TO  140  

130  CONTINUE 

140  DO  150  L=I, POINT  

IFCOKL)  .LT.THR)  GO  TO  160 

15C  CONTINUE  . . ... 

160  L=<INT<FLOAT(L)/2.0))*2 

THIS  SECTION  DIVIDES  THE  DIFFERENCE  BETWEEN  THE  CHOSEN 
BEGINNING  AND  END  POINTS  OF  THE  DENSITIES  BY  30  TO  GET 
THE  LOOP  INCREMENT  INDEX,  I,  SO  AS  TO  CALCULATE  30  ROC 
POINTS  OVER  THE  SIGNIFICANT  PORTIONS  OF  THE  DENSITIES 

I=INT(FL0AT<L-K)/3D.*C.5)  — _ _ 

IF(I.EQ.O)  1=1 

PRINT*, "ROC  PARAMETERS  -,K,L,I  _ _ „ _ 

PRINT*,-  - 

N=C  _ _ . _ _ _ 

DO  200  H:K,L,I 

N=N»1  . . 

PFA (N)=3PAL (00, M, POINT) 

PO(N)=GRAL  (01, M, POINT)  _ . 

200  CONTINUE 

PFA(N*1) =PO(N*l) =0.0  . . 

PFA(N*2)=0 .25 

PO(N*2)  = 0. 16667  . _ . 

• PLOT  THE  RECEIVER  OPERATION  CURVE  (ROC)  I PD  VS  PFA 

* 

CALL  LINE(PFA,P0,N,1,5,J)  _ 

300  CONTINUE 

CALL  PLOTE(N) 

STOP 

END  . _ 
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FAST  FOURIER  TRANSFORM  SUBROUTINE  . . _ 

SUBROUTINE  FFT (A,M,N,IS1 

A = COMPLEX  INPUT  ARRAY  OF  SAMPLES 
M s LOG2(N) 

N = 2 M =NUMBER  OF  SAMPLES 
IS  = -1»  FORWARD  TRANSFORM 

= ♦1,  INVERSE  TRANSFORM  (UNNORMALIZEO) 

COMPLEX  A»U,W,T 
DIMENSION  A(N) 

DATA  PI/ 3. 14159265/  . 

NV2=N/2 

NM1=N-1  ..  - 

Jsl 

DO  30  1=1, NHl  . _ _ 

IFCI.GE.J)  GO  TO  10 

T=A(J)  _ . 

A(U)=A(I) 

A<I)=T  

10  K=NV2 

20  IF(K.GE. J)  GO  TO  30 
J=J-K 

K=K/2  . 

GO  TO  20 

30  J=JtK  

DO  80  L=lfM 

LE=2**L  _ 

LEl=LE/2 
U=CMPLX( 1.0,0. » 

IF<IS)  40,50,50 

40  M=CMPLX«C0S(PI/LE1»  ,-SIN(PI/LEin  

GO  TO  60 

50  W=CMPLX(COS<PI/LEl),SIN(PI/LEl)> 

60  DO  80  J=1,LE1 

DO  70  I=J,N,LE  ...  

IP=I*LE1 

T*A(IP)»U  

A(IP)=A( I) -T 

70  A<I)=A(I»»T  . . 

80  U=U*H 

RETURN  _ 

END 
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Appendix  D 

The  following  program  listing  is  typical  of  the  FORTRAN 
programs  used  to  generate  receiver  performance  curves  for  the 
linear  approximate  processor.  Subroutines  not  listed  here 
are  listed  in  Appendices  B aind  C . This  program  was  also 
used  to  generate  receiver  performance  curves  for  the  two  ad- 
hoc  linear  processors  after  modification  of  the  processor 
output  density  function  calculation  subroutine. 


( 
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PROG'»a*^  GRAPH  (IMPUTjO'JTPUT,  PLOT) 

INTFGFR  POINT 

DIHFMSION  Y(60  2),  00(G0  2),ni(Fi3?) 

DIMENSION  PFA  (40)  ,!»n  (U3) 

3ATA  FAHINfOMIN," A0,00/0.0f 0.0, 0.25,0. 166S7/ 

* 

• POINT  = NUMBER  OF  DATA  POINTS  IN  EACH  DENSITY 

• INDEX  J = NUMBER  0^  DATA  SETS 

• 

READ* , POINT 

CALL  PLOT(0.  ,-3.0  ,-3) 

CALL  PLOT(0.,2.5, -3) 

CALL  PLOT(0.0, 6.0,3) 

CALL  PLOT(4.0,6.0 ,2) 

CALL  PLOT(4.0,0.0,2) 

CALL  PLOT(4.0, 6.0,3) 

CALL  PLOT(0.0, 0.0 ,2) 

CALL  AXIS(0.,0  .,26HPR0'’ABILITY  OF  FALSE  AL  ARM, -26 , 4. 0 , 

♦ 0.0, FAMIN', PAD) 

CALL  AXIS(0.,0  .,24MPP0BABILITY  0="  DETECTION  , 24 , 6 . 0 , 90  . , 
♦DMIN,nO) 

CALL  SYMBOL(0.25,6.4,0.14,24MR'^CEIVER  OPERATING  CURVE, 

♦ 0.0,24) 

DO  3G0  J=l,3 
READ* ,ETAO,ETA1,RHO 

* 

* COMPUTE  HO  AND  HI  DENSITIES 

* 

010  = 1.0-  (1.0 /P HO)  ♦((ETAl/RHO  )-pTA0)/2.0 
OIl=(RHO-1.0)*(ETAl-(R'HO*p‘^fln))  /2,Q 
PRIN-^*,"  OIO=  ’SDIO,”  011=  ",0I1 
PRIN*^*,**  •• 

IF  (DIO. LT. 0.0)  DIQ=APS(OIO) 

IF  (Dll. LT. 0.0)  DI1=ABS(0I1) 

00  100  1=1, POINT 
Y{I)=FL0AT(I)-1. 

DO (I)=XCHI (Y(I) ,ETA0, DIO) 

D1(I)=XCHI  (Y(I),ETA1,DI1) 

100  CONTINUE 

CALL  SCALE(00,4.0 ,POINT,l) 

CALL  SCALE(01,4.0 ,oOINT,l) 

THP=nO(POINT^2) 

« 

• FIND  THE  STARTING  POINT  OF  HO  DENSITY , INDEX  K 

« 

DO  110  1=1, POINT 
IF(On (I)  .6T.THR)  GO  TO  120 
110  CONTINUE 

120  Ks<INT(FLOAT(I)/2.0) )*2 
IF(K.LE.O)  K=1 
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♦ FIND  END  POINT  OF  HI  DENSITY,  INDEX  L 

* 

THR=ni(POINT+2) 

DO  ITO  1=1, POINT 
IF(D1 (I) .GT.THR)  GO  TO  140 
130  CONTINUE 
140  DO  irO  L=I, POINT 

IF(D1 (L)  .LT.THR)  GO  TO  160 
150  CONTINUE 

160  L=(INT(FL0AT(L)/2.C) )*2 

* THIS  SECTION  DIVIDES  THE  DIC-f^RENCE  BETWEEN  THE  CHOSEN 

* BEGINNING  AND  ENT  JOINTS  OF  THE  DENSITIES  BY  30  TO  GET 

* THE  LOOP  INCREMENT  iNDf^X,  I,  SO  AS  TO  CALCULATE  30  ROC 

* POINTS  OVER  THE  SIGNIFICANT  PORTIONS  OF  THE  DENSITIES 

I=INT  (FLOAT  (L-K)/ 30.  ♦•0.5) 

IFCI.EO.O)  1=1 

PRINTS, “ROC  PARAMETERS  ••,K,L,T 
PRINT*,  ••  •• 

N=0 

DO  200  M=<,L,I 
N=N*1 

PFA(N>=GRAL(D0 ,M, POINT) 

P0(N)=GRAL(D1,M,30INT) 

200  continue 

PFA(N*1)=PO(N+1)=0.0 
PFA(N+2)=0.25 
P0(N*2) =0.16667 

» 

* PLOT  THE  RECEIVER  OPERATION  CURVE  (ROC)  I PD  VS  PFA 

* 

CALL  LTNE(PPA,PO,N,l,5,J) 

300  CONT-^NUE 

CALL  PLOTE(N) 

STOP 

END 


THIS  FUNCTION  COMPUTES  A POINT  ON  THE  LINEAR  PROCESSOR 
OUTPUT  DENSITY 

FUNCTION  XCHI(Y,ETA,PI) 

DOUBLE  V 

V=-(r/OI»ETA/2.0) 

Z1=(1.0/ABS(OI).)*OFXP(V) 

Z2=FI0(SORT( (2.0*Y*ETA)/OI)) 

XCHI=71*Z2 

RETURN 

END 
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Appendix  E 

Equal  Receiver  Performance  Analysis 

This  appendix  demonstrates  analytically  that  the  non- 
linear approximate  processor  and  the  linear  approximate 
processor  have  equal  performance.  For  simplicity,  a single 
channel  processor  with  equal  conditional  variances  is 
assumed . 

Using  Eqs  (3^)»  (35) i and  (36)  the  performance  param- 
eters probability  of  detection  and  probability  of  false  alarm 
can  be  defined  for  the  non-linear  approximate  processor. 

These  parameters  are  given  by 
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where  K = c - e = Vin^-Vm^/a^  and  w is  the  threshold. 

Similarly,  using  Eqs  (3^)»  (35) » and  (38)  performance 
parameters  for  the  linear  approximate  processor  are  given  by 
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where  K'  = c^-e^/4  = nij^-mQ/4a^  and  w'  is  the  threshold. 
If  the  variable  transformation  given  by 


z 


4(VmJ-Vin^ 


(59) 


is  applied  to  the  expression  for  Pj)  for  the  linear  approx- 
imate processor  given  by  Eq  (57)  one  obtains 
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= s 
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(60) 


where  t is  the  transformed  threshold;  t = ^(Vmi-Vmo)^^^-^^ 
Comparing  this  expression  with  the  non-linear  approximate 
processor  expression  for  given  by  Eq  (55) > it  is  observed 
that  the  two  integrals  are  of  the  saume  form  and  are  equal 
when  w = t = 2( Vnii -Vino')Vw ' /(^ai -mQ )' • 

Applying  the  transformation  given  by  Eq  (59)  to  the 
expression  for  Pp^  for  the  linear  approximate  processor 
yields 
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where  t is  the  transformed  threshold  defined  above.  This 
expression  is  of  the  same  form  as  the  non-linear  approximate 
processor  expression  for  Pp^  given  by  Eq  (56).  These  expres- 
sions for  Pp^  are  equal  v;hen  w = t,  the  identical  result 


F 


obtained  above  for  probability  of  detection. 

Since  there  is  a one-to-one  correspondence  between  w* 
and  t,  corresponding  thresholds  would  generate  identical  re- 
ceiver operating  points  for  both  the  non-linear  approximate 
processor  and  the  linear  approximate  processor.  Because 
individual  operating  points  are  identical,  identical  receiver 
operating  curves  would  be  generated  for  both  processors  over 
the  entire  range  of  thresholds.  This  verifies  the  results  of 
the  nvimerical  analysis  presented  in  Chapter  IV.  While  not 
proven  here,  the  numerical  results  indicate  that  equal  re- 
ceiver performance  between  processors  can  also  be  expected 
for  the  more  general  case  of  unequal  conditional  variances, 
although  the  necessary  variable  transformation  is  much  more 
complex  than  the  transformation  presented  here. 
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